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ABSTRACT 



The shape of the spectral energy distribution of active galaxies in the EUV- 
soft X-ray band (13.6 eV to 1 keV) is uncertain because obscuration by dust 
and gas can hamper our view of the continuum. To investigate the shape of the 
spectral energy distribution in this energy band, we have generated a set of pho- 
toionization models which reproduce the small dispersion found in correlations 
between high-ionization mid-infrared emission lines in a sample of hard X-ray 
selected AGN. Our calculations show that a broken power-law continuum model 
is sufficient to reproduce the [Ne V]i4.32^m/[Ne III], [Ne V]24.32/im/[0 IV]25.89/xm 
and [O IV]25.89//m/[Ne III] ratios, and does not require the addition of a "big 
bump" EUV model component. We constrain the EUV-soft X-ray slope, ccj, to 
be between 1.5 - 2.0 and derive a best fit of ~ 1.9 for Seyfert 1 galaxies, 
consistent with previous studies of intermediate redshift quasars. If we assume a 
blue bump model, most sources in our sample have derived temperatures between 
Tbb — 10^-^^ K to 10^-^ K, suggesting that the peak of this component spans a 
large range of energies extending from ~ A600 A to A 1900 A. In this case, the 
best fitting peak energy that matches the mid-infrared line ratios of Seyfert 1 
galaxies occurs between ~ A700-A1000 A. Despite the fact that our results do 
not rule out the presence of an EUV bump, we conclude that our power-law 
model produces enough photons with energies > 4 Ry to generate the observed 
amount of mid-infrared emission in our sample of BAT AGN. 

Subject headings: AGN: general - galaxies: Seyfert - X-rays - IR 
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Introduction 



Active galactic nuclei (AGN) produce an extraordinary amount of energy which is 
thought to be primarily ther mal emission from an optically thick accretion disk fueling a 
supermassive black hole (e.g. 



Peterson et al. 



2004j ). Accretion can produce a characteristic 



feature in th e ultraviolet continuum o f most AGN known as the "big blue bump" (bbb, 



Shields 



1978 



Malkan fc Sargent 



around the Lyman limit (e.g., 



19821) and it has been argued that its energy peak lies 



Zheng et al. 



1997 



Teller et al. 



2002 



Shang et al. 



20051). The 



problem is that in AGN the soft X-ray-extreme-ultraviolet spectral energy distribution 
(SED) can be hidden because of dust and gas, thus making a direct view of the full energy 
distribution impossible. A common approach is to study high redshift quasars where the 
rest-frame extreme-ultraviolet (EUV, A ~ 100-912 A) can be observed from the ground. 
However, for the local universe at z < 0.05 we need to rely on space-borne telescopes and/or 
other indirect methods to uncover the shape of the SED in this particular range of energies. 



Within the paradigm of the unified model of AGN (e.g., lAntonuccil Il993l ) both Type 1 
and Type II AGN are intrinsically the same, with their observed difference being due to 
the visibility of the central engine. An optically-thick toroidal structure surrounds the 
central source, located at distances just beyond the broad line region (BLR). Emission 
that arises further out in the narrow line region (NLR), hundreds of parsecs in extent, is 
expected to be fairly unobstructed and isotropic, i.e., independent of our viewing angle. 
However, absorption, dust extinction and star formation may hamper t he use of soft X-r ay, 



optical and infrared obser vables on these larger-than-torus scales (e.g. 



1994 



Kraemer et al. 



20101). 



Mulchaey et al. 



One way to bypass the above limitations is to focus on high-ionization (E > 40 eV), 
mid-infrared emission lines. To be produced, they require photons with energies between 
~40-100eV (EUV-Soft X-ray), and so we can assume that they are primarily associated 
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with the AGN, in other words, they are less contaminated by star formation processes. 
Also, they are less sensitive to dust extinction than optical lines. In this regard, it is 



possible to fit these observed lines vi a photoionization mode 



continuum in this energy band (e.g.. 



Alexander et al. 



1999 



s in order to constrain the AGN 



20001 ). Among the mid-infrared 



properties of Seyfert galaxies, the correlations between high-ionization emission lines, 
[Ne III] A15.56 fim, [O IV] A25.89 fim and [Ne V] AA14.32/24.32 ^m, are strikingly tight 



with a dispersion of less than ~ 0.2 dex, both in flux and luminosity 



2007 



2010 



Dudik et al 



2009 



Gou 



Pereira-Santaella et al 



ding fc Alexander 



2009 



Melendez et al. 



(see 



2008b 



Goriian et al. 



Weaver et al. 



20101 ). The [Ne V] emission, given the high energy require to 



be p roduce (E >98 . 9 eV) , has been unambiguou s 



see 



Satyapal et al. 



2007 



Abel fc Satyapal 



2008 



y identified with the presence of an AG N 



Bernard- S alas et al. 



2009 



Weaver et al. 



2010l ). therefore, the strong correlation between this emission and the lines produced from 



(E > 54.9 eV 



Melendez et al. 



sugges t that these lines 
fl2008al ) probed the utility 



transitions of [Ne III] (E > 41.0 eV) and [O IV 
are produced primarily by the AGN. Moreover, 
of high-ionization mid-infrared lines as an indicator of AGN power. They found a tight 
correlation, with a dispersion of ~ 0.3 dex, in Seyfert 1 galaxies, between the [O IV] and 
the X-ray 14-195 keV continuum luminosities. A weaker correlation was found for Seyfert 
2 galaxies, which we believe is due to the effect of Compton Scattering in the 14-195 keV 



(e.g., 



in some Sey 


fert 2 


Riebv et al. 


2009; 



galaxies. This result 


las since been confirmed 


Diamond-Stanic et al. 


2009; 


Weaver et al. 


2010) 



Our aim is to recover the hidden portion of the spectral energy distribution of AGN 
between the EUV and the soft X-ray (13.6 eV to 1 keV). To do this we have performed 
extensive photoionization calculations to characterize the small dispersion found in the 
correlations between high-ionization, mid-infrared emission-line ratios. We compare our 
calculations with emission line ratios from Spitzer/IRS spectra of the Swift Burst Alert 
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Telescope (BAT) AGN samplqj of the local universe, presented in 



Weaver et al 



torn . All 



of our calculations have been carried- o ut by the photoionization code CLOUDY, version 



08.00, last described by 



Ferland et al. 



fll998h 



2. The Models 



We have generated two sets of photoionization models to study the i ntrinsic SEP of 



thejonizing continuum. The first assumes a broken power law, as us ed by 



Melendez et al. 



1998 



(|2008af) and similar to that suggested for NGC 5548 and NGC 4151 (IKraemer et al. 
2000), of the form oc u'", with a= 0.5 below 13.6 eV, at from 13.6 eV to 1 keV and 0.8 
at higher energies. To examine specifically the EUV/soft X-ray band, we varied at from 1.0 
to 2.5. Our second model is a parametrization of the AGN continuum that combines an 
X-ray power law and a "blue bump" component, in the form: 



U = u'^- exp {-hu/kTBB) exp {-kTin/hu) + ai/"% (1) 

where auv is the low-energy slope of the Big Blue Bump continuum, Tbb is the UV 
bump cut-off temperature and kTj^ is the infrared cut-off temperature of the big bump, 
corresponding to 9.1 fim. For the X-ray power-law, ax is the slope of the X-ray component 



and t he coefficient a is adjusted to produce the correct UV to X-ray slope (ITananbaum et al. 



19791 ). aox, as defined by 

/.(2keV) 



/.(2500A) 



403.3"°^ (2) 



^The Swift/BAT covers the whole sky at (l-3)xl0~^^ergs cm~^ s~^ and represents a 
complete sample of active galaxies including Compton-thin AGN that were missed from 
previous X-ray surveys in the 2-10 keV band because of their high column densities {Nh ~ 
10^4 cm-2 ). 
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For the "big bump" model we ado pted a value for the UV to X-ray spectral slope of 



a. 



-1.4 



Zamorani et al. 



AGN fiLusso et al 
to auv = -0.5 (e.g. 



19811 ). in good agreement with recent studies of Type I 



2010r). We fixed the low-energy slope of the big blue bump continuum 



Francis et al. 



1993h . For the X-ray 



of Or 



-1. 



Korista et al. 



his parametrization is similar to that of 



power-law, we adopted a value 



Mathews fc Ferlandl fll987l ) and 



(119971 ). Finally, for this model we considered five values for the UV bump 



cutoff temperature, from Tbb = 10^'^^ K to 10®" K. Figure [T] compares the SED for our 
different photoionization models. 

To investigate the physical conditions in the emission line regions for the [O IV], [Ne III] 
and [Ne V] lines, we started with a simple, single-zone, constant density model. The logs of 
abundances of Oxygen and Neo n relative to H by number are -3.31 and -4.0, respectively. 



where the O abundance is from 



AUende Prieto et al. 



(120011 ) and the Ne abundance is from 
Holwegerl (l200ll ). We considered hydrogen column densities of lO^^cm"^ and lO^^cm"^. 



Despite the fact that dust is likely mixed in with the emission-line gas in the NLRs of 



Seyfert galaxies (e.g.. 



Kraemer fc Harrington il986l : INetzer fc Laor 



19931 ). the tightness of 



the mid-infrared correlation suggest that extinction in t he mid-infrared, on average, is not 



an important bias in our sample of X-ray selected AGN (jWeaver et al. 



2010|). 



We generated a grid of photoionization models varying the total hydrogen number 
density , (rih). and the ionization pa rameter, U, where the ionization parameter U is defined 



as see 



Osterbrock fc Ferland 



20061): 



U 



1 



hu AttR^chh 



(3) 



AirR^cnH 

where R is the distance to the cloud, c is the speed of light and Q{H) is the flux of ionizing 
photons. In the next section we will compare these models with the line ratios observed 
from the BAT AGN sample. 
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3. Comparison with Observations 



3.1. Density in the Emission Line Regions 



The use of emission hnes that are close in ionization potential but different critical densi- 



ties c an be used as a density diagnostic of photoionized plasmas (e.g., 



Osterbrock fc Ferland 



20061 ). In particular, the mid- infrared [Ne V] ratio , [Ne V]i4.3?um /fNe V]9A.?.9.jim, has been 



Dudik et al. 



2007 



Tommasin et al. 



20101 ) 



used density diagnostic in active galaxies (e.g. 
with the caveat that a wavelength-dependent mid-infrared extinction can affect this ratio, 
as suggested by the nur nber of Seyfert galaxie s that have [Ne V] ratios below the theoretic al 



low d ensity limit (e.g.. 



and 



Tommasin et al 



2010 



Weaver et al 



2010 



20101). However, t he m id-infrared extinction curve derived 



Li &: Draine 



Pereira-Santaella et al. 



Weingartner fc Draind (l200l[ ) 



200ll ) fails to predict the required 



ain the observed [N e V] ratios 



20101 ). Therefore, 



Li fc Draind ( 120011 ) (e.g, see Figure 16 in 
dust extinction at Ai4.32/im/NH and A24.32/im/NH to exp^ 
that fall below the low density limit in the BAT sample ( IWeaver et al. 
since the relation between the [Ne V] ratio and extinction is not fully understood, we 
minimized the possible effect of reddening toward the NLR by focusing on the [Ne V] 
ratio observed in Seyfert 1 galaxies as detailed multi-wavelength studies on Seyfert galaxies 
sugge st more extinction of emission from the NL R in Seyfert 2 galaxies than in Seyfert Is 



(e.g. 



Melendez et al. 



2008a; 



Kraemer et al. 



properties of the BAT sample of AGN (jWeaver et al. 



2011 



. Followin g the work on the mid-infrared 



20101 ) we found an average neon ratio 



observed in Seyfert 1 galaxies of [Ne V]i4.32^m/[Ne V]24.32/im = 0.95 ± 0.28 (see Table [T]). 

We used the optimizqj command in CLOUDY to vary Uh and U to reproduce the 
observed neon ratio. From this, and by using our "power-law" model with = 1.5 for a 



^For a full description on the optimization methods see CLOUDY documentation. Hazy 1, 
at http://www.nublado.org/ 
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total hydrogen column density of uh = lO^^cm"^, we obtained a value of Uh = 10^'^^ (cm~^^ 
and logf/ = —1.37. For our "blue bump" model, assuming a cutoff tempe rature of 
Trr = 10^'^ K such that the energy of the UV bump peaked at ~ 22 eV (e.g. 



Korista et al. 



19971 ). we obtained a value of Uh = 10^'^^ (cm~^) and logf/ = —1.51. For a total hydrogen 
column density of Nh = 10^^ cm~^ we obtained a value of rih = 10^'^^ (cm~^) and 
logU = —1.31 for our power-law models and, from our "blue bump" model, assuming 
a cutoff temperature oi Tbb = 10^'^ K, we obtained a value of rih = 10^'^^ (cm~^) and 
logf/ = —1.30. In Figure [2] we present the neon ratio versus the total hydrogen number 
density for these models for a range of ionization parameters and EUV/Soft X-ray slopes. 
Figure [2] shows that, close to the observed values, the neon ratio is not strongly sensitive to 
the SED of the ionizing radiation. 

From the different models and continuum approximations we obtained a total hydrogen 
number density of ~ 10'^ (cm~^). One must note that this value represents the density 
only of the [Ne V] emission-line gas. However, this density is in agreement with the 



va 
in 



ues found from the pho toionization modeling of the [O IV] emitting region presented 



Melendez et al. 



( l2008al ). suggesting that the [Ne V] and [O IV] emitting regions are 
essentially co-located . On the other hand, this density is smaller than that found by 
Baskin fc Laorl ( l2005l ) for the [O III] emitting region in their single-zone approximation, i.e.. 
Tie ~ 10^-^^ (cm-3). However, they found that the observed [O HI] A5007, [O HI] A4363 and 
H/3 lines can also be fitted using two-zone model comprised of a a dense {rie ~ 10^ cm~^) 
inner zone, dominated by the [O HI] A4363 emission, and a low density, more extended zone 
(rig ~ 10^ cm~^), where most of the [O HI] A5007 arises. Our density value is similar to that 
used in their low-density zone, although the densities in their two-zone approximation are 
based on fixed assumed values, whereas our density is derived from the fit to the observed 
[Ne V] ratio. 
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3.2. The EUV-Soft X-ray slope 

Following the results from the previous section, and in order to keep the number of free 
parameters to a minimum, we used the observed [Ne V] ratios to fix the hydrogen density 
in our photoionization models to Uh = 10^ (cm~^). Next, in order to maximize each of the 
mid-infrared emission lines, we proceeded to find the range in ionization parameter where 
the predicted ionic column densities for the Ne III, Ne V and O IV peak (see Figure [3]). 
Results show that the ionization parameter falls in the range —2.0 < logU < —1.5 for 
Nh = lO^^cm"^, while it is in a higher ionization range, logf/ ~ —1.0, for Nh = lO^^cm"^. 
One must note that the peak of the ionic column density depends on the SED (see Figure [5]). 
Also, because we are interested in the shape of the ionizing continuum, we expect to achieve 
the clearest constraints from radiation-bounded models q The set of radiation-bounded 
models was obtained by initially assuming a hydrogen column density of Nh = lO^^cm"^. In 
Figure H] we compared the total Hydrogen density as a function of ionization parameter and 
UV slope for our radiation bounded, "power-law" models. It is clear that a column density 
of Nh = lO^^cm"^ is enough to generate an optically thick, radiation-bounded cloud for the 
set of ionization parameters and hydrogen density derived from the previous considerations, 
except for = 1.0 where we had to extend the column density to Nh = lO^^cm"^. Using 
this set of ionization parameters, and = 10^ cm~'^, we compared our model predictions 
with the observed [Ne V]u.32fxm/[^e III], [Ne V]24.32/im/[0 IV] and [O IV] /[Ne III] ratios 
from the BAT sample. Overall, these results suggest that the [O IV] and [Ne V] come from 
a hig her ionization state and lower density regions than that suggested by the optical [O III] 



lines (IBaskin fc Laoi 



2005 



Melendez et al 



2008a). 



^Specifically, in a radiation-bounded model, all of the incident radiation has been absorbed 
inside the cloud, i.e., the cloud is optically thick to the ionizing radiation. 
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3.3. The Power-Law model 

As we mentioned before, high-ionization mid-infrared ratios can be used to investigate 
different excitation mechanisms, thus different shapes of the ionizing continuum. They 
provide information about the relative flux density of ionizing photons with energies 
between ~40-100 eV. In Figure |6] we compared the observed [Ne V]i4.32;im/[Ne 111] and 
[Ne V]24.32/im/[0 IV] ratios found in the BAT sample (see Table [1]) against those predicted 
by our photoionization calculations. Because of the proximity of wavelengths and equal 
aperture extraction for the emission lines in each ratiqj, this comparison is unaffected 
by dust extinction and aperture effects. Figure |6] shows two sets of "power-law" models, 
corresponding to Nh = 10^^ cm^^ and Nh = 10^^ cm~^, for a range of the EUV/soft X-ray 
from = 1.5 — 2.0. We show that, for a column density of Nh = 10^^ cm^^, the cloud is 
optically thin to the ionizing radiation, as the models underpredict the amount of [Ne III] 
emission. Specifically, this result describes a matter-bounded component, in which some 
fraction of the ionizing radiation passes through the gas unabsorbed. Overall, a softer 
ionizing continuum will result in a lower [Ne V]/[Ne 111] ratio, because there are fewer 
photons capable of ionizing Ne^^. 

Figure |6] shows that, independent of the ionization parameter, most of the observed 
ratios are matched by models with EUV/Soft X-ray spectral indices of 1.5 < «j < 2.0. 
Using the 95% confidence interval for the actual mean for the [Ne V]i4.32/im/[Ne III] 
and [Ne V]24.32/im/[0 IV] ratios observed in Seyfert 1 galaxies in our sample, we found 
a smaller range of spectral indices, i.e., 1.7 < < 2.0, with ~ 1.9 providing the best 
fit for the median observed values. A histogram of the EUV/Soft X-ray spectral indices. 



^The [Ne V] i4.32yum and [Ne III] fiuxes were extracted from the Spitzer/YR^ in the short- 
high order (SH, 4.7"x 11.3"). On the other hand, the [Ne V]24.32/im and [O IV] emission 
were extracted from the Spitzer/\Yi^ in the long-high order (LH, ll.l"x 22.3") 
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derived empirically from the position of each source in Figure |6l is presented in Figure [71 
From their analysis of Hubble Space Telescope jHS T) Faint Obje c t Spe ctrograph (FOS) 
observations of a sample of 101 quasars at z > 0.33, IZheng et al.l (Il997l ) found that the 
power-law index in the EUV between 350 and 1050 A is a = 1.96 ± 0.15. They also found 
that, for a sub-sample of 60 radio-loud sources, the EUV spectral index is a ~ 2.2, while 
for a sub-sample of 41 radio-quiet it is a ~ 1.8. As shown in Figure El most of the BAT 
AGN can be fitted with an EUV/Soft X-ray spectral index of a < 2.0, consistent with the 
radio-quiet nature of most Seyfert galaxies. 



From their analysis of HST spectra of a sample o 



184 quasars at z > 0.33 



Telfer et al. 



Zheng et al. 



(Il997l ). This sample 



(120021 ) found a slightly harder EUV continuum than 
included observations with the HST FOS, Goddard High Resolution Spectrograph and 
available Space Telescope Imaging Spectrograph (STIS) data. For their complete sample, 
they found a spectral index of a = 1.76 ± 0.12 between 500 and 1200 A . For a subsample 
of radio-loud objects, they found a = 1.96 ± 0.12, while for a radio-quiet subset, they found 
a = 1.57 ± 0.17. These values are roughly consistent with our range of EUV/Soft X-ray 
indices and the spectral index associated with the median values for the Seyfert 1 galaxies 
found in our sample, a ~ 1.9. Again, this supports the idea that a single power-law can 
represent the EUV/soft X-ray continuum in our sample of BAT AGN. One must note 
that star formation contamination in the [Ne HI] emission can bias the EUV/Soft X-ray 
slope towards a steeper power-law index. In other words, star formation contamination 
in the [Ne HI] emission will result in a lower [Ne V]i4.32/im/[Ne HI] ratio, which could be 
wrongly associated with a softer EUV continuum. However, the tight correlation found in 
the [Ne III]-BAT relationship suggests that, on average, there is no strong enhancement 
due to star formation in the [Ne HI] emission in our sample, although we cannot rule out 
relatively minor effects in individual objects. Despite the fact that our spectral index has 
been determined over a relatively small range of energy (~40-100 eV), these results support 
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the validity of our method, given the agreement with values derived from different spectral 
analyses, covering a wider range of energies. 

On the oth er hand, ana l ysis o f the Far Ultraviolet Spectroscopic Explorer {FUSE) 



archival data by 



Scott et al. 



( I2OO4I ) revels a significantly harder composite spectral energy 



distribution, a = 0.56lo;28; than the spectral indices discussed above and no evidence of a 



UV spectral index break. One must note that, in their sample, some individual sources do 
show a spectral break in the far-ultraviolet (A ~ 910-2000 A). It should be noted that this 
is an heterogeneous sample of AGN, which were observed with FUSE because they were 
known to be bright in the UV. The authors reported a correlation between the spectral 
slope and AGN luminosity; in other words, low- luminosity AGN tend to show a harder 
ionizing continua. They suggested that the harder spectral index they derived is due to the 
relatively large fraction of low luminosity AGN in the Ff/S'i? samplqj. In contrast, we found 
no evidence of such flatter indices in our sample. Specifically, in Figure |8] we compare our 
derived values for the spectral index and the 14-195 keV luminosities, from which it can 



be seen that there is no depend ence between the 



the BAT sample. Furthermore, 



Shang et al. 



lUV/Soft X-ray slope and luminosity in 



(120051 ) found that most objects in their FUSE 



sample exhibit a spectral break arou nd 1100 A , in agreement with previous iJ>ST results 



(Zheng et al. 


1997; 


Teller et al. 


2002) 



2OO2I ). Their analysis is based on a heterogeneous sample 



of 17 low-redshift AGN (z< 0.5), with quasi-simultaneous spectrophotometry including 



FUSE, HST and Kitt Peak National Ob servatory o 
from 900 - 9000 A. Given the fact that Scott et al. 



j servat ions covering a rest wavelength 



( 12004] ) looked for possible systematic 



errors that could bias their results, the source of this discrepancy is unclear. While it is 
possible that the spectral break is luminosity dependent, another possible explanation is 



^The median value for their FUSE sample is log ALuoo = 45.0 ('ergs/s) versus log A Lnnn = 
45.9 (ergs/s) for the i^ST sample of active galaxies presented by iTelfer et al.l (120021 ) 
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that intrinsic reddening, if not fully correct ed for, may lead to an incorrect determination 



of the spectral break and continuum shape (IShang et al. 



20051). 



Figure [9] shows the predicted and observed [Ne V]i4.32/im/[Ne III] and [O IV]/[Ne III] 
ratios. Contrary to the ratios presented in Figure [6l the [O IV]/[Ne III] ratio could be 
susceptible to reddening and aperture effects. In this regard, these effects will have different 
manifestations in the observed [O IV] /[Ne III] ratios. For example, the [O IV] /[Ne III] 
ratio will have larger values if extinction is stronger for the shorter-wavelength line. On the 
other hand, aperture effects will tend to decrease the observed ratio. As seen in Figure [9l 
when using the [O IV]/[Ne III] ratio, we recovered the same range for the EUV/soft X-ray 
slope as those derived from the [Ne V] 14.32^111/ [Ne III] and [Ne V]24.32/im/[0 IV] ratios 
(see Figure |6]). In other words, the [O IV] /[Ne III] ratios are fit by the range of spectral 
indices found from the other ratios, suggesting that reddening and aperture effects are not 
dominant (or, perhaps, cancel e ach other out) in the BAT sample of AGN, in agreement 



with the discussion presented in 



Weaver et al. 



(120 101 ). 



3.4. The AGN Parametrization 

Figures [lUHII] show predictions for the [Ne V]i4.32^m/[Ne III], [Ne V]24.32/im/[0 IV] 
and [O IV] /[Ne III] ratios from a set of models generated with the parametrization of the 
AGN continuum including a UV bump, with cutoff temperatures from Tbb = 10^'^^ K to 
10®" K, and column densities of = 10^^ cm~^ and = 10^^ cm~^. Once again, at 
Nh = lO^^cm"^ the cloud is optically thin to the ionizing radiation, thus underpredicts 
the [Ne III] emission. It is clear that the AGN parametrization with a UV "bump" 
temperature higher than 10^'^ K fails to predict most of the observed ratios, even in our 
radiation-bounded models. We find that the emission-line ratios for a majority of sources 
can be fit by models with UV bump temperatures of Tbb = 10^'^^ K to 10^'^ K, or. 
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alternatively, peaks at ~ 600 1900 A. The best fit for the 95% confidence interval for 

the actual mean occurs between ~ 700 1000 A. These results are in agreement with the 

FUV break observed in the compos i te spectra of AGN , in the range of 1000-1100 A (e.g., 



Zheng et al. 



1997 



Telfer et al. 



2002 



Shang et al. 



20051). 



This range of temperature is slightly lower than the typical ter nperatures used to 



characterize an average quasars continuum when modeling the BLR (IKorista et al. 



19971) 



However, the chosen form for the parametrization of the AGN continuum presented in 



Korista et al. 



(119971 ) focused on luminous quasars while our sample is dominated by lower 
luminosity sources. If the ionizing continuum is from an accretion disk, one would expect 
more luminous objects to be hotter for the same black hole mass, and, therefore, the lower 
range of temperature found in our sample is not unexpected. However, these results may 
also argue in favor of a different, intrinsic SED photoionizing the BLR. This intrinsic SED 
may be characterize by a bbb with a higher "blue bump" temperature, e.g., Tbb > 10^'^ K, 
than that found from our modeling of the NLR emission. In this interpretation, material 
located between the BLR and the NLR may absorb the intrinsic UV component. However, 
from our calculations, we do not find any clear evidence of a modification of the SED by an 



intervening medium. 



4. Conclusions 

In order to investigate the shape of the spectral energy distribution in the range 13.6 eV 
to 1 keV, we have generated photoionization models to characterize the small dispersion 
found in correlations between high-ionization mid-infrared emission-line ratios in AGN 
This is the band where dust extinction and the line of sight absorber can make it virtually 
impossible to have a direct view of the continuum generated by the nuclear engine. Our set 
of radiation bounded models show that a broken power-law continuum model is sufficient 
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to reproduce the observed mid- infrared [Ne V]i4.32/im/[Ne III], [Ne V]24.32/im/[0 IV]25.89/im 
and [O IV]25.89/im/[Ne III] ratios. From this we constrained the EUV/Soft X-ray spectral 
index, a^, to be between 1.5 - 2.0, with ~ 1.9 corresponding to the median value obse rved 
galaxies. Our result ag rees with previous studies (e.g.. 



in Seyfert 



Telfer et al. 



2002 



Shang et al 



Zheng et al. 



1997 



20051 ) and holds for the wide range of BAT luminosities and 
type I and II objects in the local universe, whereas previous studies of SED tend to focus 
on UV selected objects at a different redshift range. Despite the fact that our spectral 
index has been estimated for a relatively small range of energy (~40-100 eV), these results 
probe the validity of our method to uncover the AGN continuum in this energy band. From 
the comparison between the photoionization models and observations we derived empirical 
values for the EUV/Soft X-ray slope for each source and found no dependence between the 
EUV/Soft X-ray slope and the 14-195 keV luminosity in our sample. 

We also find that our power-law model can produce enough photons at energies >4 Ry 
to generate the amount of mid-infrared line emission observed in the BAT AGN sample. 
However, this result does not rule out the presence of a bbb. When we do add a "big 
bump" in the UV to an X-ray power-law, the majority of sources required "blue bump" 
temperatures of Tbb = 10^'^^ K to 10^ '' K, suggesting that the peak of this component 
lies in a wide range of energies that can extend from ~ 600A to 1900 A. The best fit 



agreement with i: 


revious estimates { 


Telfer et al. 


2002 


Shang et al. 


2005) 



~ 700-1000 


A, in 


Zhene et al. 


1997; 



20051 ). This range in temperatures is lower than the typical 



UV b ump peak assumed to characterize an average quasar continuum when modeling the 



BLR ( iKorista et al. 



19971 ). However, if the ionizing continuum is from an accretion disk, 
one expects more luminous objects to be hotter for the same black hole mass. Therefore, 
the lower range of temperature found in our sample is not unexpected. 
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Photon Energy (eV) 

Fig. 1. — Comparison between the different continuum spectral energy distribution con- 
sidered in tfie present work. For tfie sake of comparison, for the "blue bump" model we 
present two values for the UV bump cutoff temperature 10^'^° K and 10^ ° K. Again, for 
comparison purpose, wc used a "power-law" model with a value for the EUV/Soft X-ray 
index of ctj = 1.5. The three continua shown are normalized to have the same integrated 
luminosity between 1 Ryd< hu < 7.354 x 10^ Ryd, however, because we are focusing on the 
modeling of emission line ratios, instead of of absolute line fluxes, the chosen normalization 
is arbitrary. 
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Fig. 2. — Upper panel. The [Ne V] ratio versus the hydrogen number density, uh for 
the "power-law" model with a, = 1.5 for different ionization parameters, logU — —0.5 
and logU = —3.0. Lower panel. The [Ne V] ratio versus the hydrogen number density, 
uh, for the "power-law" model for different spectral slopes with a fix ionization parameter, 
logC/ = —1.5. The solid horizontal line represents the observed average ratio and the 1-a 
{dashed lines) found for Seyfert 1 galaxies. 
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Fig. 3. — Predicted ionic column densities for Ne III {circles), Ne V (triangles) and O IV 
(squares), for two set of models, as a function of ionization parameter U, and two values of 
the column density, lO^^cm"^ and lO^^cm"^, left and right panels, respectively. The sohd 
line represents the "power-law" model with = 1.5 and the dashed line represents the "big 
bump" parametrization with the peak of the UV bump at « 44 eV, or alternatively a UV 
bump cutoff temperature of 10^ K. 




Fig. 4. — Comparison for the radiation bounded models between the total Hydrogen col- 
umn density and ionization parameter for the "power-law" model for different values of the 
EUV/Soft X-ray index, a^. One must note from this comparison that, within the range 
of ionization parameters that reproduce the observed ratios, our radiation bounded mod- 
els yield to Hydrogen col umn densities higher than that found in typical NLR conditions. 



logA^H = 21 (cm-2) (e.g., 



Kraemer et al 



20001 ) 
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Fig. 5. — A comparison between the ionization parameter associated with the peak of the 
ionic column density versus the EUV/Soft X-ray spectral slopes for the "power-law" model. 
From this comparison is clear that the peak on the ionic column density depends on the 
shape of the ionizing continuum. 
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Log [NeV] /[OIV]^ 



Fig. 6. — Comparison between the predicted and observed [Ne V]i4.32;im/[Ne III] and 
[Ne V]24.32/im/[0 IV] ratios. Seyfert 1 galaxies are presented as black circles, Seyfert 2 galax- 
ies are red triangles, blue squares represent the newly detected BAT AGN and green stars are 
LINERs. We present two set of "power-law" models corresponding to log Nh = 21 (cm^^) 
{solid lines) and log Nh = 22 (cm~^) [dashed lines) for a range of = 1.5 — 2.0. For the 
set of models presented in this comparison we show a range of ionization parameters of 
—2.5 < logU < —1.5, in steps of 0.5 dex and 0.1 dex for log Nh = 21 and log Nh = 22, 
respectively. All these models corresponds to a hydrogen number density of Uh = 10^ (cm^). 
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Fig. 7. — Distribution of EUV/Soft X-ray slopes, ai, derived from our analysis of the entire 
BAT sample. One must note that the spectral index was derived empirically from the 
position of each source in Figure [6l therefore, each slope has an aj±l error. For any given 
source located between two consecutive power-law models, e.g., different spectral slope, we 
selected the flatter index as the empirical estimate. 
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Fig. 8. — Comparison between the spectral index, ai, and the 14-195 keV luminosities 
where it can be seen that there is no dependence between the EUV/Soft X-ray slope and the 
14-195 keV luminosity in our sample. One must note that the spectral index was derived 
empirically from the position of each source in Figure |6l therefore, each slope has an ctjil 
error. For any given source located between two consecutive power-law models, e.g., different 
spectral slope, we selected the flatter index as the empirical estimate. 
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Fig. 9. — Comparison between the predicted and observed [Ne V]i4.32;im/[Ne III] and 
[O IV] /[Ne III] ratios. Seyfert 1 galaxies are presented as black circles, Seyfert 2 galaxies 
are red triangles, blue squares represent the newly detected BAT AGN and green stars are 
LINERs. We present two set of "power-law" models corresponding to log Nh = 21 (cm^^) 
{solid lines) and log Nh = 22 (cm~^) [dashed lines) for a range of = 1.5 — 2.5. For the 
set of models presented in this comparison we show a range of ionization parameters of 
-2.5 < logU < -1.0, in steps of 0.5 dex and 0.1 dex for log Nh = 21 and log Nh = 22, 
respectively. All these models corresponds to a hydrogen number density of Uh = 10^ (cm^). 
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Fig. 10. — Comparison between the predicted and observed [Ne V]i4.32/im/[Ne III] and 
[Ne V]24.32/xm/[0 IV] ratios. Seyfert 1 galaxies are presented as black circles, Seyfert 2 
galaxies are red triangles, blue squares represent the newly detected BAT AGN and green 
stars are LINERs. We present two set of "blue bump" models corresponding to logA''^:^ = 
21 (cm^^) (so/zd lines) and log Nh = 22 (cm~^) {dashed lines) for five values of the UV bump 
cutoff temperature between Tbb = 10^'^^ K and 10^ ° K. For the set of models presented 
in this comparison we show a range of ionization parameters of —2.5 < logf/ < —1.0, in 
steps of 0.5 dex for both log Nh = 21 and logNfj = 22. All these models corresponds to a 
hydrogen number density of Uh = 10^ (cm^). 
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Fig. 11. — Comparison between the predicted and observed [Ne V]i4.32/im/[Ne III] and 
[O IV] /[Ne III] ratios. Seyfert 1 galaxies are presented as black circles, Seyfert 2 galaxies 
are red triangles, blue squares represent the newly detected BAT AGN and green stars are 
LINERs. For the sake of simplicity, we present two set of "blue bump" models corresponding 
to logA'^i^ = 21 (cm^^) {solid lines) and log Nh = 22 (cm~^) {dashed lines) only for three 
values of the UV bump cutoff temperature, Tbb = 10^'^^ K, 10^ '^° K and 10^ ° K. For the 
set of models presented in this comparison we show a range of ionization parameters of 
-2.5 <\ogU < -1.0, in steps of 0.5 dex for both log Nh = 21 and logA^^ = 22. All these 
models corresponds to a hydrogen number density of Uh = 10^ (cm^). 
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Table 1. Statistical Analysis Between Seyfert 1 and Seyfert 2 Galaxies 



Seyfert 1 Seyfert 2 





Measurements 




Standard 


Measurements 




Standard 


Pk-S 




Available 


Mean 


Deviation 


Available 


Mean 


Deviation 


(%) 


[Ne V]i4.32/im/[Ne V]i4.32//m 


29 


0.95 


0.28 


30 


1.00 


0.38 


87.8 


[Ne V]i4.32/im/[Ne III] 


32 


0.47 


0.22 


33 


0.41 


0.20 


79.8 


[Ne V] 24. 32^111 /[O IV] 


31 


0.34 


0.19 


30 


0.28 


0.06 


57.0 


[O IV] /[Ne III] 


38 


1.58 


0.87 


33 


1.65 


0.77 


91.4 



Note. — The last column, Pk-Si represents the Kolmogorov-Smirnov (K-S) test null probability. Upper limits for the 
[Ne V] fluxes are not included. This table also includes information about the numbers of Seyfert 1 and Seyfert 2 galaxies, 
mean values and standard deviations of the mean for the measured quantities. For the [Ne V]i4.32;um/[Ne III] ratios the 

95% confidence interval for the actual mean is 0.3904 — 0.5465 for Seyfert 1 galaxies, and 0.3417 thru 0.4831 for Seyfert 2's. 
For the [No V]24.32^m/[0 IV] ratios the 95% confidence interval for the actual mean is 0.2717 thru 0.4089 for Seyfert 1 
galaxies and 0.2469 thru 0.2964 for Seyfert 2's. Finally, for the [O IV]/[Ne III] ratios the 95% confidence interval for the 
actual mean is 1.295 thru 1.867 for Seyfert 1 galaxies, and 1.375 thru 1.924 or Seyfert 2's. 



